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Abstract. A general formula is obtained from which the Madelung type constant: 
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extensively used in the finite-size scaling theory is computed analytically for some 
particular cases of the parameters d and v. By adjusting these parameters one 
can obtain different physical situations corresponding to different geometries and 
magnitudes of the interparticle interaction. 
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In the analytic investigation of the finite-size scaling theory of systems undergoing phase 
transition the Madelung type constant 
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where I G T, d and T(a,x) is the incomplete gamma function, plays a central role. 
By adjusting the parameters d and u, one can obtain constants, describing different 
physical situations. These situations correspond to different geometries (hypercube, 
slab geometry and many others) and interparticle interaction in the system (short as 
well long-range). The particular case C(4|l), corresponding to the short-range forces, 
has been widely used in the asymptotic analysis of the finite-size properties of the 
O (n) -symmetric <£> 4 model using renormalization group treatment of static 0, |3], |J as 
well dynamic 0, |^, 0] critical phenomena. The constant C (d\v) for the long-range case 
has been obtained in the asymptotic analysis of finite-size effects of the spherical model 
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of Berlin and Kac JT], §, |], [10], |TTJ as well as the quantum tp A model [|12j in the large 
n limit. The same constant is obtained in the renormalization group treatment of the 
finite-size scaling in O (n) -symmetric systems [ O . 

The constant C(4|l) is evaluated numerically with a very good accuracy. It is found 
to be (see e.g. [fj) 

C(4|l) = -1.7650848012... tt = -5.545177444.... (2) 

It is the aim of this letter to find a general formula for the analytic evaluation of 
the constant C{d\v) and subsequently to deduce useful expression in some particular 
cases of the parameters d and v. Let us note that the integral C(d\v) has the remarkable 
symmetry property 



u v C {d\u) = ii d,2 - v C[d 



which relates the values of C (d\u) for v > | with those for v < 4. Equation (|3|) is 
obtained as a consequence of using the Jacobi identity for the sum in the integrand. 

Our key finding is that the Madelung type constant can be expressed in terms of 
the analytic continuation, over v < |, of 

c(d\v) = J- 2 »r(v)J2' r 2v , v>\, (4) 

i 

where I e Z d and the primed summation indicates that the term corresponding to I ^ 
is excluded. For some particular values of the dimension of the lattice, the d-iold sum 
can be expressed as a product of simple sums such as Ditrichlet series ||14j| . 

To show that the Madelung type constant (|1]) is related with the sum given in 
equation (01) we start from the generalized d- dimensional Jacobi identity: 



(5) 



Following reference [T5|], we multiply both sides of (|5[) by u d l 2 v 1 and integrate over u. 
Whence, we get the key identity (valid for v ^ 0, |) 

I - Z7 ^ i d - 2!/ ^ l 2u W 

* Z I 

The right hand side C(d\v) of identity (0) is a constant of integration independent 
of u. Consequently the right hand side should also be u independent. By adjusting 
the parameter u one obtains different expressions for the constant C(d\v). All these 
expressions are equivalent in the sense that they give the same 'numerical' value for 
fixed d and v. In particular, we find it is useful to obtain simple expressions for the 
constant Cid\v) corresponding to the limiting cases u —>■ oo and u — > 0. With the aid of 
the asymptotic behavior of the incomplete gamma function [I5| 

^ e -*[l-Sr± + 0{£)], X>1, 

T(a,x) = { (7) 
T(a) - $e-" [l + ^ + O (x 2 )] , x « 1. 
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it is possible to perform the evaluation of the right hand side of (|6|) and one gets ([I]) 
(valid for < v < |) in the limit u — > 0, and (f|) (valid for v > |) in the limit u — > oo. 
It is not difficult to see that the results (0) and (Q) are valid in two different intervals 
and so they complement each other. 

On the other hand the sum in the right hand side of equation (^) can be expressed 
in terms of the Epstein zeta function [|T7| 



Z 



(d } v) = J^r 2 \ lez d , u>- } (8) 
I 

which can be regarded as the generalized ci-dimensional analog of the Riemann zeta 
function ((v). In the case under consideration the Epstein zeta function has a simple 
pole at v = | and may be analytically continued in the interval v < ~. Note that, from 
the functional equation for the Epstein function |I7| , one can check easily that C (<i|z/), 
defined in (|), obeys the symmetry property (|3]). 

Using the results of JR], [15], for the Epstein zeta function (||), we get simple 
expressions, for C(<i|z/), for certain values of d: 

(a) For the simplest one dimensional case, d = 1, we obtain 

C{l\v) = 27r 1 / 2 -^i>)C(2^) ^0,1 (9) 

As a particular case, we give here the value of the constant C(l 1 1/4) = 2r(l/4)£(l/2) = 
—10.589351..., corresponding to the short-range case with v = |. 

(b) In the two dimensional case, d = 2, we obtain 

C(2\u) = 4 7 r 1 - 2i T(z,)C(z/) / W v?0,l, (10) 
where (3{v) is the analytic continuation of the Dirichlet series: 

oo 

^) = E(- 1 ) / ( 2/ + 1 ^ v >0 - 

1=0 

Note here that in the particular case corresponding to a long-range potential with v = ~, 
we get C(2|l/2) = 4 v /?C(l/2)/3(l/2) = -6.913039577.... 

(c) The constant C(d\v) in the four dimensional case, d = 4, turns out to be 

C(4\u) = 8 (1 - 4 1 -") Tr^r^K^ - l)C(iz) i/ ^ 0, 2. (11) 

Particularly we are interested in the value of the constant in the case of short-range 
interaction corresponding to v = 2. We find it to be exactly |T5 



C(4|l) = -8 In 2. 

As we mentioned above this constant has been widely used in the analytic investigation 
of finite-size scaling and its relation to numerical analysis. 

Now we turn our attention to the three dimensional case. To our knowledge there 
are no analytic expressions for it up to now. This interesting case has been investigated 
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numerically in reference [ I7| . In order to evaluate numerically the constant C(d\v), here 



we propose the following general formula 

(12) 
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obtained from equation (y) for u small enough. Equation fll2|) is a generalization of 



a three dimensional result obtained in reference [p!q] . Note that equation ([12]) is valid 



for arbitrary d and z/ 7^ 0, |, which makes it suitable for numerical evaluations of the 
constant C(d\u), especially in the cases of 3d and 5d used in finite-size scaling. 



The nature of the error term equation ( |T2| ) is such that this formula is useful 
for numerical evaluations even when the parameter u is not too small. Since C{d\v) 
with concrete values of the parameters d and v are related with finite-size properties 
of confined systems with different geometries and different types of interparticle 
interaction, we shall present some most useful values obtained numerically from equation 

(0)- 

In the case of short-range interaction two constants were used in the literature 
corresponding to different situations. For a three dimensional system confined to a 
fully finite geometry we find the value C(3| 1/2) = —8.91362917... in perfect agreement 



with references [|], ||, [18], [0|. For a system with a slab geometry, we get C(3|l) = 
—5.028978843..., coinciding with the value of reference 0. Another case of interest is 
that of system with long-range case v — 3/4 and a cubic geometry. For this case we 
have C(3|3/4) = —5.9098415587.... At the end, we quote a result for the 5c? case, it is 
C(5|l) = —4.228709895... in agreement with that of reference ||], where this constant 
was used in the investigation of a the finite scaling scaling in a five dimensional system 
confined to a cubic geometry. 

Numerical values of C(Av\v)/T{y) for v = |, |, | and the exact value for v = 1 



were recently reported in reference ||20|| . 

This work is supported by The Bulgarian Science Foundation under grant F608/96. 

Note added in proof 

After the completion of this paper we learned that in a particular form, identity |] was 
first used by Epstein in 1903. We are thankful to R.M. Ziff, who brings to our attention 



this fact and the review article on ideal Bose-Einstein gas in reference ||21|| , where the 
same type of functions appear in the case of finite volume. 
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